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ABSTRACT 



The purpose of. this study is to investigate the 
influence of non-environinental and small scale environ- 
mental acoustic factors on the variation of transmission 
loss within areas of the ocean which exhibit homogenie ty of 
large scale environmental factors. The method employed was 
to identify the influencing and non-influencing factors, 
pool over the non-influencing factors if they existed, and 
investigate the natvure and magnitude of the effect of the 
influencing factors on the variation of transmission loss. 

The results generally indicated that source/receiver 
depth, the mixed layer, geographic location, and transmit- 
ting range and frequency significantly influenced the 
variation of transmission loss. 
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I. INTRODUCTION 



In recent years, with the advent of more sophisticated 
SONAR systems, there has been an increasing emphasis placed 
on the study of the acoustic properties of the ocean. To 
be able to predict the capabilities of a particular SONAR 
system, whether it be passive or active in nature, requires 
a complete understanding of the effect of the local acous- 
tic environment on the transmission of sound through the 
water. Both theoretical and observational techniques have 
been utilized in an effort to determine the environmental 
and non-environmental factors which effect the transmission 
of sound and then to characterize this transmission as a 
function of those factors. 

One study of this nature was the AMOS project performed 
under the sponsorship of the U. S. Navy in the early 1950 's 
(ref. 5). In this study, the behavior of sound transmission 
was measured along with local environmental characteristics 
over a wide range of environmental and non-environmental 
factors. On the basis of these observations, transmission 
loss prediction equations were developed as a function of 
certain factors such as source and receiver depths, range of 
transmission, mixed layer depth, etc. However, these equa- 
tions are quite insensitive to small scale environmental 
characteristics such as inhomogenie ties in the thermal 
structure and irregularities at the air-sea and sea-earth 
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interfaces. As a consequence, transmission losses observed 
operationally are often quite variable with respect to the 
predicted losses on the basis of these equations. 

This study is an extension of the work performed in the 
AMOS experiment in that it investigates the nature of the 
small scale environmental factors which influence the 
observed temporal and spacial variability in acoustic pro- 
pogation. The characterization of this variability as a 
function of these factors will be useful in a variety of 
applications. Some examples are: 

1) Establishing some design specifications 

for new SONAR systems. 

2) Establishing SONAR detection probabilities 

for shipboard use. 

3) Simulation of environmental conditions for 

use with detection simulation models. 
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II. A MODEL OF THE ACOUSTIC ENVIRONMENT 



Within sufficiently small areas of the ocean, ’ the large 
scale environmental factors (i.e., thermal, density and 
salinity profiles) may he considered spacially consistent. 
Such areas may also be considered temporally consistent for 
short durations of time. Since general acoustic behavior 
is a fvinction of these large scale . environmental conditions, 
such an area may be considered acoustically homogeneous. 

This, however, does not imply that the acoustic behavior is 
independent of small scale temporal and spacial inhomogeni- 
eties within this large scale environmental structure. 

Hence, the environment may be characterized as a generally 
homogeneous medium which dictates the large scale acoustic 
behavior, but one having small temporal and spacial inhomo - 
genieties which cause variations in acoustic behavior. 

On the basis of this characterization, sound transmission 
loss may be modeled as a stationary stochastic process 
whose mean is a fianction of the large scale environmental 
factors, and whose autocovariance function is a function of 
the small scale environmental factors. Denote this process, 
where X(t) describes the transmission loss, by X(t);t>0. 
It is convenient to characterize this process by its mean 

/c = E(X(t)) 
and by its autocovariance function 
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C( A t) = E((X(t) - /U) (X(t + A t) - /c )) . 

It is then possible to study the nature of the variability 
which occurs in transmission loss due to small scale 
environmental aspects by considering the environmental and 
non-environmental factors which effect C( A t) . 

Two particular parameters which characterize this 
variability are C(0), the variance of the process and hence 
the amount of variability that would be experienced within 
an acoustically homogeneous region for a given source and 
receiver depths, and At*, the minimum value of At for 
which C(At) = 0, which represents the shortest time inter- 
val between independent observations. These two parameters, 
C(0) and At*, may be determined for a given process if 
that process is observed either continually or at discrete 
increments of time. Then by observing such processes under 
varying acoustic conditions (such as thermal structure, 
geographic location, time of day and year, etc.) it should 
be possible to determine how C (O) and At* vary as a 
function of the acoustic environment. 

However, it is not necessary to observe the entire 
process to determine C (O) and At*. Specifically, if 
the process of transmission loss is observed at two points 
in time separated by A t and these pairs of points are 
observed many times at varying values of At for the same 
process, or equivalently from processes with the same 



10 



autocovariance functions, then this information may be used 
to estimate C ( A t) and hence, C(0) and At*. This may 
be accomplished in the following way: 

Let 

Y(At) = X(t) - X(t+ At) 

represent the difference between two measurements in trans- 
mission loss taken in the same location but separated in 
time by At. 

Then 

E(Y( A t)) = 0 

and 

Var(Y( At)) = E(Y( A t)2) 

= Var (X(t)) + Var (X(t + At)) - 
2C (At) 

= 2 (C(0) - C( A t)) 

The value of Var(Y( At)) may then be determined for a 
given At from the observations of Y ( A t) corresponding 
to pairs of observations of transmission loss with the same 
autocovariance functions, C( A t) . C (O) may then be 
represented by 

C(0) = Var(Y( A t*))/2 
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where Var (Y( A t*) ) is the value of Var (Y( A t) ) when 
C( A t) approaches zero or when Var(Y( At)) no longer 
increases with At. At* is then the time interval where 
this first occurs. 

If pairs of observations of transmission loss cannot be 
obtained from the same process, then it is necessary to pool 
information from different processes with the same auto- 
covariance function. The problem here arises in specifying 
the range of conditions under which this may be done. The 
approach taken in this study was to characterize the local 
acoustic environment by its thermal structure where it was 
hypothesized that environments with similar thermal struc- 
tures will produce processes with similar if not identical 
autocovariance functions. This assumes that other variables 
such as transmitter and receiver depths and transmission 
ranges will remain constant. 
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III. NATURE OF THE DATA 



The data utilized in this report was made available by 
the U. S. Underwater Sound Laboratory, New London, 
Connecticut, and reflects the raw data collected during 
the Acoustic, Meteorological and Oceanographic Survey (AMOS) 
conducted from June, 1949, through April, 1953. There were 
nine cruises staged during these four years which covered 
the North Atlantic, the Norwegian Sea, and to a lesser 
extent, the Mediterranean Sea. 

Two ships were employed during each cruise where one 
acted as the transmitting platform and the other as the 
receiving platform. Each of the cruises was divided into 
several widely spaced stations which served as the focal 
points for data measurement and collection. Within each 
station, the acoustic data was collected at several trans- 
mission ranges between 800 and 30,000 yards. At each 
transmission range the transmission loss measurements were 
observed for four transmitting frequencies at various 
source/receiver depth combinations. The transducer depths 
varied from 20 to 500 feet. The specific data recorded was 
transmission loss, in decibels (db) , as a function of 
source/receiver depths and range between ships. Table X 
lists the information that was recorded. 

Accuracy of the data becomes particularly important in 
consideration of the accuracy of any results obtained 
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through the use of that data. The following list gives an 
approximate measure of the accuracy of the AMOS data: 



1) time - to the nearest minute 

2) range - to the nearest 5 yards 

3) transducer depths - to the nearest foot 

4) water depth - to the nearest 10 fathoms 

5) latitude and longitude - to the nearest minute 

although it is possible that the two ships 
may have drifted as much as a mile or two 
from the recorded position 

6) BT pattern code - reflects the gross thermal 

conditions closest to the time that the 
data was recorded. 
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IV. ASSUMPTIONS 



As indicated in the introduction, this paper deals with 
determining the variation of transmission loss by examining 
the variation of the change in transmission loss as a 
function of the elapsed time between a pair of observations. 
Ideally it would be best if this change in decibel loss 
could be measured from a pair of observations with identi- 
cal source and identical receiver depths, and with the same 
horizontal transmitting range, but separated by an interval 
of time, A t. Although data of this nature was not availal^le 
from the AMOS experiment, pairs of observations of db loss 
were obtained having the following characteristics: 

1) One of the observed pair measured the 

transmission loss of a signal generated 
at depth d^^ and received at depth d 2 
while the other measured transmission loss 
of a signal generated at depth d 2 and 
received at depth d]^. 

2) Identical transmitting ranges between the 

source and receiving ships was not always 
available within a pair of observations. 

Figures (1) and (2) depict the nature of an ideal pair 
of observations and that of the available data. 

The fact that the change in transmission loss was 
measured using transmission paths that 'crossed' instead of 
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Figure 1. Ideal Transmission Paths 



Position 




being 'parallel' was considered not to be a serious depart- 
ure from the ideal conditions (refs. 3 and 5). This follows 
from the assumption of a generally constant horizontal 
thermal environment, and therefore an acoustically 
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homogeneous environment over the range of operation of the 
two ships on a given station, and theoretical considera- 
tions which postulate’ that transmission loss over the same 
path from different directions will be the same. It is 
foreseen, however, that this approximation to the ideal case 
of parallel transmission paths will increase the observed 
■Var(Y(At)) for small At due to a contribution caused 
by small scale spacial variability, and may also cause a 
slight underestimation of At. However, the estimation 
of C(0), the measure of total variation experienced within 
a generally homogeneous acoustic environment, will not be 
effected. 

The second data characteristic of range differences 
between the two observations in a pair was treated in the 
following manner. The observed change in transmission loss 
was adjusted for the effect of this range differential on 
the basis of a piecewise linear approximation of transmis- 
sion loss vs. horizontal transmission .range. This linear 
approximation seemed reasonable in the interest of keeping 
the range adjustment computation simple and yet maintaining 
sufficient accuracy. Plotting transmission loss vs. trans- 
mission range for each source/receiver depth combination at 
each station revealed that the resulting relationships 
exhibited the same general shape and that these relation- 
ships could be approximated by two linear segments with the 
break point at 5,000 yards. 
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As the slopes of these two linear segments are a function 
of the local environmental conditions, and thus varied from 
station to station, it was necessary to estimate the slope 
for each segment at each station for each source/receiver 
depth coiribination. This was done on the basis of the trans- 
mission loss observed at the various ranges occupied on a 
given station. Then the change in transmission loss 
observed at a particular station and source/receiver depth 
was adjusted using the corresponding slope of transmission 
loss vs. transmission range in the following manner: 

^^^adjusted ~ ^ observed “ ^ ^ ^ range 

where A db represents the change in transmission loss, b 
is the slope of the corresponding linear segment, and 
A range is the range differential within a pair of obser- 
vations. 

In some cases insufficient data was available to 
estimate the slope of one or both segments at a given 
station and source/receiver depth combination. In such 
cases the corresponding adjustments were made on the basis 
of a slope which was typical of other stations on the same 
and other cruises. Table I depicts a representative cross- 
section of the average slopes from the available data. 

Listed below is a summary of the typical slope values selec- 
ted for use within this study where sufficient data was not 
available for slope calculations. 
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Range Interval 


Slope Value 


Applicable 

Cruises 


Above 


.0015 


5, 7, 8, 9, 10, 


5,000 yards 




11, 12 


Below 


.0045 


7, 8, 9, 10, 12 


5,000 yards 


.0100 


11 




.0280 


5 



Admittedly this adjustment procedure is subject to 
error. In particular, the error incurred by using a linear 
approximation in the first place, and secondly in the esti- 
mation of the slope of the transmission loss vs. transmission 
range relationship. It was felt, however, that in general 
the total error inherent in this procedure could be 
restricted to at most i.5db by making no adjustments with 
range differentials within pairs of observations of more 
than 500 yards. Errors of such magnitude were considered 
acceptable in light of the original accuracy in measuring 
transmission loss in the AMOS experiment. In actual 
practice, this error was further limited by the fact that 
few adjustments in change in transmission loss were made for 
range differentials greater than 250 yards. 
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V. RESULTS 



A. METHODOLOGY OF INVESTIGATION 

The purpose of this research is to investigate the 
nature and magnitude of the variability in sound propogation 
loss as characterized by C( A t) by investigating the 
variability in the change in transmission loss as a function 
of time between observations. The calculation of trans- 
mission loss as a function of time, denoted by Y( A t) , 
was obtained from the AMOS data in all instances where a 
pair of transmission loss observations were taken under the 
following conditions; 

1) The compared observations were taken at 

the same station. 

2) The range differential between observations 

was no greater than 500 yards. 

3) The source depth of one transmission equalled 

the receiver depth of the other, and 
vice-versa, thus resulting in a compari- 
son of observations with similar but not 
identical transmission paths. 

These calculations of change in transmission loss over 
common paths were indexed by the following environmental 
and non-environmental factors; 

1) Source and receiver depths 

2) Transmission range 
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3) Transmitting frequency 

4) Water depth 

5) Mixed layer depth 

6) Geographical location 

7) Time interval between corresponding observations 

8) Time of day and year 

The investigation of variability in sound transmission 
was confined to the influence of the above factors. To 
determine this influence using the available data, it was 
necessary to pool data collected under varying environ- 
mental conditions. Such pooling is only applicable, however, 
over a range of environmental conditions where the acoustic 
variability is characterized by a common autocovariance 
function. Hence, it is necessary to determine how to iden- 
tify the range of conditions which lead to a common auto- 
covariance function. In particular, it is necessary to 
determine which factors, both environmental and non- 
environmental, influence the nature and magnitude of the 
variation of propogation loss in the ocean and then pool 
over the entire range of those factors which did not so 
influence this variation. 

The procedure followed was to plot Var(Y( At)) vs. 

At using the entire sample set of data and thus represent- 
ing the entire range of environmental and non-environmental 
conditions encountered. This entire sample of data was 
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then subdivided into sub-samples on the basis of the 
source/receiver depth combination and Var(Y(.A t) ) vs. 

At was plotted for each sub-sample. By this procedure, 
it was determined that the nature and magnitude of 
Var(Y( At)) depended upon the sub-sample (i.e., source/ 
receiver depth combination) considered. Hence, it was 
deemed inappropriate to pool the data over source/receiver 
depth combinations when considering the influence of other 
factors on Var(Y( A t) ) . 

To determine if the thermal structure influenced 
Var ( y ( At)), each of the above sub-samples was further 
subdivided on the basis of the location of the receiver anci 
source depths with respect to mixed layer depth. On the 
basis of this subdivision, it was determined that the posi- 
tioning of the source and the receiver with respect to the 
mixed layer depth also influenced the nature and magnitude 
of Var (Y ( A t) ) . 

The influence of geography and transmission range were 
also investigated in the same manner by subdividing the 
sample data on the basis of these factors. As it was 
necessary to keep the respective sub-samples as large as 
possible when considering a given factor, it was found that 
certain sub-samples should not be further subdivided. Such 
was the case for some of the less representative source/ 
receiver depths. 
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By proceeding in this manner, it was determined that, of 
the factors being considered, those which were judged to 
have a significant influence on the nature and magnitude of 
Var(Y( At)) were: 1) source/receiver depth, 2) location 

of the source and receiver with respect to the mixed layer, 
3) geographic location, and 4) transmission range. It was 
also determined that the transmitting frequency had a marked 
effect on the magnitude of Var (Y ( At)). The actual nature 
of these influences will be discussed later. 

It should be emphasized that Var(Y( At)) is twice 
the value of C(0) and thus Var(Y( t) ) , the variation of 
transmission loss. 

B. LIMITATIONS EFFECTING RESULTS 

There are two considerations which tended to limit the 
results of this paper. The first was that the amount of 
data available to estimate Var(Y( At)) as a function of 
At was limited in some of the finer subsamples defined in 
the previous section. Hence it was found necessary, in 
order to obtain a reasonably precise estimate of 
Var(Y( At)), to pool the available data into several con- 
secutive 15 minute intervals. 

The second consideration was that, due to physical 
limitations imposed on observation of the data during the 
AMOS experiment for a given soiarce/receiver depth combina- 
tion, the majority of the available data was restricted to 
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at most three 15 minute intervals. These three intervals, 
however, were not necessarily the same for each S 9 urce/ 
receiver depth combination, but tended to the smarll values 
of At for closely spaced source and receiver depths and 
to the large values for widely spaced source and receiver 
depths . 

The implication of these data limitations was that A t 
could not be determined with a precision greater than 15 
minutes o However, this was thought not to be a serious 
limitation. A possibly more serious limitation is that, 
because of the dispersion of the data as a function of At, 
the value of Var (Y ( At)) cannot be estimated with equal 
precision from interval-to-interval and from depth 
conibination-to-depth combination. In particular, for a 
given source/receiver depth combination, there may be little 
or no data available for the interval when Var (Y ( A t) ) 
reaches its maximum value. In general it is felt that this 
did not happen, although there is no way of verifying this 
at the present time. 

C. FINDINGS 

1. Distribution of Transmission Loss 

As suggested in the introduction, knowledge of the 
nature of the variability in transmission loss may be 
utilized to produce a realistic simulation of the acoustic 
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environment. To do this it is necessary to know not only 
the magnitude of this variation, but also its distribution. 

The typical assumption that this distribution is normal 
is evaluated on the basis of the available data for a 
representative cross-section of siibsamples of the data. 

The basis for this examination is the Kolomogorov-Srairnov 
one-sample goodness-of-f it test (ref. 2) . Table II gives 
the results for all subsamples considered, and the range of 
conditions which define the tested subsample. Also given 
are the maximum absolute difference between sample 

cumulative distribution and a normal cumulative distribution 
with mean zero and sample variance of the subsample. The 
critical value has the same connotation usually associated 
with hypothesis testing and represents the level of signifi- 
cance at which the null hypothesis of no departure from 
normality would be rejected. 

The general conclusion on the basis of these results is 
that the normal distribution is an adequate approximation 
for the distribution of transmission loss. 

2. Factors Effecting the Variation of Transmission Loss 

a. Source/Receiver Depth 

The first attempt at determining which environ- 
mental and non-environmental factors influenced the nature 
and magnitude of Var (Y ( At)) was in the area of 
source/receiver depths. The available data was subdivided 
on the basis of several source/receiver depth combinations. 
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The results thus obtained proved to be quite noisy and it 
was decided to employ the same procedure for cruise 10 only 
since the thermal structure appeared more consistent on that 
cruise. The results of this procedure seemed satisfactory 
for exposition purposes and are summarized in Table III. 

It is interesting to note that not only does Var(Y(At)) 
increase as time increases, but that this variation also 
increases with an increasing depth differential between 
source and receiver. Also noted was that Var(Y( A t) ) 
tended to decrease as the source depth increased. The 
latter is probably attributable to surface effects, 
b. Mixed Layer Depth 

Once it had been decided that the different 
source/receiver depth combinations each influenced 
Var(Y( At)) differently, the next step was to consider 
the position of the source and receiver with respect to the 
depth of the mixed layer. This was accomplished by further 
subdividing the data into two subsamples on the basis of 
whether both the source and receiver were in the mixed 
layer, or at least one was below the mixed layer. The 
rational behind this subdivision was to determine if the 
variability of transmission loss was different when the 
transmission path was above the thermocline or across the 
thermocline. 

Table IV shows the results of this investiga- 
tion and clearly indicates that Var(Y(At)) is 
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considerably less within the mixed layer. The differences 
observed may in part be a fionction the existence of internal 
waves . 

c. Geographic Location 

Geographic location was foiond to have a 
significant effect on the variation of transmission loss. 

To facilitate investigating the nature of this effect, areas 
were selected which were thought to have approximately the 
same large scale environmental characteristics. The areas 
chosen were: 



1) 


50-80 


N, 


all longitudes 


2) 


0-50 


N, 


0-60 W 


3) 


0-50 


N, 


60-90 W 



Four source/receiver depth combinations were then 
picked which would ensure that sufficient pairs of observa- 
tions would be available within each of these depth combina- 
tions. The siibsamples have been further subdivided accord-? 
ing to the position of the source and receiver with respect 
to the mixed layer. Table V depicts Var(Y( A t) ) for the 
geographic areas and the four depth combinations selected. 
The variation of transmission loss for each combination 
pooled over all latitudes and longitudes has also been 
included in the table for comparative purposes. 

As indicated by the entries in Table V, there 
does not appear to be any consistently predictable direc- 
tion or magnitude associated with Var(Y( At)) across the 
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three geographic locations. However, the variations do seem 
to change significantly both from the "all latitude- 
longitude" case, and across geographic areas. This would 
seem to substantiate the hypothesis that geographic loca- 
tion is a significant factor effecting Var(Y( At)). The 
data in Table VI, which breaks the variation of trans- 
mission loss down into 15 minute intervals, also appears 
to substantiate this contention. 

d. Transmitting Frequency 

The effect of transmitting frequency on 
Var(Y( At)) is apparent; increasing the frequency in- 
creased the magnitude of the variation. Two frequencies 
were selected to investigate this hypothesis, 8 and 25 KC. 
Table VII is a summary of Var (Y ( A t) ) for both frequen- 
cies for several source/receiver depth-geographic location- 
mixed layer depth combinations. The increase in Var(Y(At)) 
is evident. 

e. Transmission Range 

The transmission range appeared to have a 
significant effect on Var(Y( At)). Table VIII shows this 
variation for range intervals of 0-10 Kyds and 10-30 Kyds 
for pairs of observations at the 8 KC transmitting frequency 
pooled appropriately according to source/receiver depth, 
geographic location and mixed layer depth. An examination 
of this data suggests that Var (Y ( At)) was generally 
lower for greater ranges and that estimates of variation 
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in the mixed layer was less consistent than for those 
transmissions crossing the layer. 

3. Investigation of Large Values of Change in 
Transmission Loss 

Because of the limited amoxant of data in some 
subsamples considered, the estimator of Var(Y( At)) was 
quite sensitive to a not uncommon occurrence of dispropor- 
tionately large values of the change in transmission loss. 

In an attempt to determine the nature of the circumstances 
which lead to such observations, a sample of observed 
changes in transmission loss whose absolute value exceeded 
9db were examined as to the conditions under which they were 
observed. This sample consisted of all instances for 
20/100, 50/100, 50/250 and 100/500 foot source/receiver 
depths at 8 KC transmitting frequency, and 20/100 and 100/500 
foot source/receiver depths at 25 KC transmitting frequency. 
Table IX is a summary of these data. 

On examination of these data it is concluded that, 
except for a marked effect of frequency and a suggested 
effect of time of day, no common single or set of circum- 
stances is related to the large changes in transmission 
loss. This does not mean that such a set of circumstances 
does not exist, but that it is not apparent what such a set 
would be based on in consideration of these data and this 
specific set of indices. 
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This latter point has lead to the consideration 
that there may be a more appropriate set of indices with 
which a measured change in transmission loss may be classi- 
fied. The nature of such a set of indices is suggested by 
considering the results of the present study. That is, that 
Var (Y ( At)) varied with the proximity of the source or the 
receiver to the surface, the depth differential between the 
source and receiver, the relationship of the source and 
receiver to the mixed layer, geographic location, etc. 

These factors are exactly those which influence the mode of 
transmission of a signal in the general transmission loss 
equations obtained in the AMOS study (ref. 5) . 

In that study it was determined that transmission 
loss depended on the mode of transmission. In particular 
on whether transmission was by 

1) Direct path 

2) Single surface bounce 

3) Multiple surface bounce 

4) Leakage by diffraction or scattering 

5) Depressed sound channel 

6) Bottom bounce 

On the basis of the results of the present study, 
it appears that the nature and magnitude of the variability 
in transmission loss might also vary with transmission mode. 
Further, it is envisioned that subdividing the available 
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data on the basis of transmission mode will yield a more 
uniform distribution of At for a given subsample and, in 
general, increase the amount of data available for the esti- 
mation of Var (Y ( At)) for any particular interval of At. 



31 



I 



VI. CONCLUSIONS AND RECOMMENDATIONS 



The aim of this study was to determine the environmental 
and non-environmental factors which influence the nature 
and magnitude of the variability in transmission loss, and 
then to characterize this variability as a function of these 
influencing factors. 

On the basis of this study it was found that the 
variation in transmission loss was influenced by the depth 
of the source and receiver of the transmitted signal, the 
relationship of the source and receiver depths with the 
mixed layer depth, geographical location, and transmission 
range. 

The nature of this influence was that the magnitude of 
variability increased with increasing transmitting frequency, 
with increasing depth differential between the source and 
receiver, with the increasing closeness of the source and 
receiver to the air-sea interface, with decreasing trans- 
mission range, and when transmitting across the thermocline 
as opposed to transmitting exclusively within the mixed 
layer. 

The time interval between uncorrelated observations was 
found to be shortest near the surface as opposed to at 
depth, to increase with increasing source/receiver depth 
differential, and to decrease when transmitting across the 
thermocline as opposed to transmitting entirely within the 



mixed layer. 
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The magnitude of the variability in transmission loss 

experienced under the varying conditions ranges over one 

2 

order of magnitude from 10 db in the near horizontal trans- 
mission in the mixed layer at 8 KC to 100 db crossing the 
thermocline in the 25 KC case. Such figures are only esti- 
mates and as such are subject to statistical variation, but 
the trends in magnitude of the variability appear to be 
consistent in the various cases considered. Therefore, the 
magnitudes of the variability determined in this study may 
be used as approximations of the true variability which will 
be ejq)erienced in the corresponding environmental situation. 

In conclusion, it is observed that the factors which 
influenced the nature and magnitude of the variability in 
transmission loss are exactly the same as those found in the 
AMOS study to influence the mean transmission loss (ref. 5) . 
In particular, the factors of source and receiver depth, 
mixed layer depth, transmission range and surface tempera- 
ture were used to determine the mode of transmission of a 
given transmitted signal. The transmission loss expected 
was then calculated using an equation peculiar to this speci- 
fic transmission mode. 

In a like manner it was suggested that the nature of the 
variability in transmission loss may be more adequately 
classified and the magnitude of the variability more pre- 
cisely specified if the observed change in transmission 
losses were further indexed and examined by transmission 
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mode. It is therefore recommended that this be the next 
area of investigation in the study of the nature and magni- 
tude of the variability of transmission loss in the ocean. 
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APPENDIX A 



I. DESCRIPTION OF DATA 
A. TAPE RECONSTRUCTION 

The information as received from the U. S. Underwater 
Sound Laboratory was in the form of 9, seven-track tapes 
containing appropriate special, temporal, and acoustic 
information. These tapes have now been consolidated into 
two nine-track tapes suitable for mounting on the standard 
IBM tape drive units. One tape contains only the bathy- 
thermograph (BT) information while the other contains all 
other necessary data including acoustic measurements 
pertinent to this report. Table X displays the information 
that is recorded on this latter tape. The following inforraa- 
tion has been omitted from each record on the revised 
acoustic tape, but is still available on the original seven- 
track tapes; 

1) Ambient noise indicator 

2) Wind force (beaufort) 

3) Unknown contents following water depth entry 

The acoustic tape contains data for cruises 5, 7, 8, 9, 
10, 11, and 12. Each record is 72 characters in length 
although each only utilizes 60 character spaces. The addi- 
tional 12 spaces have been left blank to allow room for any 
additional information deemed appropriate and to facilitate 
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placing this data on a user disc if desired. No special 
records have been added to indicate the end of a cruise or 
station. It should be noted that all field entries are in 
terms of integer numbers with the exception of fields 5, 

24, and 27 which will require character format (A-format 
under FORTRAN IV) . 

The data for cruises 2 and 3 were deemed inappropriate 
and therefore not included on the acoustic tape. The 
reason was that the source/receiver depths and two of the 
four transmitting frequencies do not correspond to the 
requirements of this study. That is, both cruises show 
source depths set at a constant 15 feet while the receiver 
depth is varied. Consequently, none of this data would be 
accepted by the computer program. The program requires the 
comparison of one record with source depth at d^^ feet and 
the receiver depth of d 2 feet, with another record having 
the source at depth d 2 and receiver at depth d^^. 

B. COMBINING ACOUSTIC DATA 

An examination of the original tapes revealed that 
either 1, 2, or 3 acoustic measurements were taken for each 
particular location, time, and range. Hence, some method 
must be employed to combine these sets of data in order that 
they be consistent. One approach would have been to take an 
average each time a grouping of two or three measxirements 
was encountered. Another approach, and the one sxibsequently 
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adopted, would be to simply sum the acoustic data for 
groupings of three measurements, sum and multiply by 3/2 
for groupings of two, and triple any single entries. The 
resulting figures would then be entered on the revised 
acoustic tape in place of the old transmitting frequency 
data, giving a single entry for each particular location, 
time, and range. This latter, method was utilized in the 
interest of conserving valuable computer time, and in 
consideration of computational and programming ease. In 
summary, all acoustic field entries at each of the four 
transmitting frequencies must be divided by 3 before they 
can be used in any computations. 

C. BATHYTHERMOGRAPH PATTERN CODE AND ITS USE 

The bathythermograph (BT) pattern code is a two digit 
number and will be found in fields 18 and 20 on the acoustic 
tape. 

Interpretation of the code has been simplified for use 
within this study. Figure (3) indicates the general nature 
of the thermal structure. To use the figure, enter with 
the first digit of the code across the top. This gives a 
picture of the approximate temperature gradient down to 100 
feet (associated depths may be read in the '1st digit' 
column) . Next, enter with the second digit. This gives a 

picture of the approximate temperature gradient from 100 to 
250 feet. Connecting the two pictures thus obtained results 
in a complete recap of the thermal conditions from 0 to 250 
feet. 
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Depths (ft) Pattern Code Digit 

1st digit 2nd digit 12345678 




Figure 3. - BT Pattern Code Interpretation 

This BT pattern code information was used to establish 
the position of the mixed layer. 
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APPENDIX B 



I. EXPLANATION OF COMPUTER PROGRAM 



This section has been included in the report in order 
that the program may be used as a base for further studies 
in the field of temporal and spacial effects on the varia- 
tion of decibel loss underwater. The entire program has 
been written and compiled under FORTRAN IV, G-level. 

A. GENERAL EXPLANATION 

The program is designed to read one data record at a 
time from the acoustic tape, check the appropriate data 
against the parameters established by the input array, and, 
if possible, pass the entire record to subroutine COMPl. 
Subroutine COMPl then separates the records passed to it 
into four categories according to range between the trans- 
mitting and receiving ships, the range break point, and 
source and receiver depths. For example, if the source/ 
receiver depths are either 20/100 or 100/20 feet and the 
range break point is 5,000 yards, the categories are: 



Category 



Soxarce/Receiver 

Depth 



Range Interval 



III 



II 



I 



20/100 

20/100 

100/20 

100/20 



- 5,000 yards 
^5,000 yards 
5- 5,000 yards 



IV 



5 , 000 yards 
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Once this data has been collected for a single station, the 
appropriate slope information is computed by subroutine 
RGRESS and the next station is then processed through the 
main program, COMPl, and RGRESS. 

At the end of each cruise, the processed records are 
sent to subroutine COMP. This routine computes the change 
in transmission (db) loss in a sequential manner. That is, 
a record can never be used in a computation more than once. 
The computation of A db loss proceeds in the following 
manner : 

1) At any given time, two adjacent records are 

' considered. 

2) The program computes the difference between 

ranges of the two records. If this difference 
is greater than the specified amount 'Z', then 
no computation is made and the next data set 
is picked up. 

3) If the range difference is less than or equal 

to the specified amount, Z, then the change 
in db loss is calculated. 

Of the. two records being considered at any particular time, 
the decibel loss of the second is adjusted in the direction 
of the first record with respect to the difference in range 
associated with the two records. This is the point at which 

the slope information passed from subroutine COMPl is util- 
ized. The only exception to this computational procedure 
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is when one range falls below the break point and the other 
above. In that event, both db losses are adjusted to the 

range break point. Once a computation has been completed, 
the change in decibel loss and the time difference between 
the two records are stored in appropriate arrays. At the 
end of each cruise, this information may be plotted at the 
users option. Also within this subroutine are two arrays 
which collect the same information for all cruises to 
facilitate a cumulative plot. 

At the end of the tape, the program calls subroutines 
SIGMA and BTLAYR. SIGMA simply sorts the A db losses with 
respect to their associated time increments into 5 and then 
15 minute intervals, computes the variance within each 
interval, and plots the variance vs. the change in time. 

It also calculates the variance and standard deviation over 
all data passed to it at the end of each tape iteration. 

BTLAYR is designed to consider all comparisons which 
lead to computation of Adb loss in subroutine COMP, and 
separate the pairs of observations into one of three cate- 
gories: 1) above the layer, 2) transmitter and receiver 

both below the layer, and 3) all transmissions with the 
receiver and source neither above or below the thermocline. 
Once all records have been categorized, the change in db 
loss vs. change in time may be plotted at the users option. 
Next subroutine SIGMA is called, which computes the 
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appropriate output as described above for each of the three 
transmission categories. 

Subroutine SIME solves simultaneous linear equations 
necessary to subroutine RGRESS. OSPLOT is a standard IBM 
plotting package and is used exclusively for plotting within 
the program. 

B. INPUT SPECIFICATIONS AND DESCRIPTION 

An input array has been established in order to lend a 
certain amount of flexibility in filtering data through 
specific parameters. This information must be punched on 
computer cards for inclusion at the end of the program and 
will be either integer or alphabetic in nature. There must 
be six entries on each card with the last letter of the 
numeral entered successively in columns 6, 12, 18, 24, 30, 
36, 42, 48, 54, 60, and 66. The last 9 entries on the last 
card may be left blank. Table XI fully explains what 
information is required and gives a general description of 
the input variables. The input data is read into the NK 
array. Note that NK(15) and NK(17) require a character 
format. Tables XII and XIII describe the arrays and vari- 
ables bsed within the computer program. 

C. READING DATA FROM THE ACOUSTIC TAPE 

The acoustic tape data format is described by Table X. 
There is one peculiarity which exists within the formatting 
structure. Note that the range between the two ships is 
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divided into two separate field entries. In this case it 
is necessary to compute the range in the following manner 
(tape records are loaded into the IN array) : 

Range = 1N(9) + 1,000 X IN (8) 

All data on the tape is integer in nature except fields 5, 

6, 24 and 27, which are alphameric. 

D. ESTABLISHING THE MIXED LAYER DEPTH 

Within each pair of compared data sets are four BT 
pattern codes, one per ship per record. Since there exists 
a change in time and a range differential between data sets, 
it seemed appropriate to combine the depths associated with 
the pattern codes in some manner. One acceptable method and 
the one adopted, takes an average of the associated thermal 
layer depths across records first for the transmitting ship, 
and then the receiving ship. This averaging is done in the 
same manner for determining both the top and the bottom of 
the thermocline. Once this has been accomplished, the mini- 
mum of the two top of the layer, and the maximum of the two 
bottom of the layer averages are selected to represent the 
upper and lower bounds of the thermocline respectively. 

As an example, consider some BT pattern code XjlX 2 » 
where it was desired to determine both the top and the 
bottom of the thermocline. Computing the top and bottom of 
the thermocline was done in the following manner: 
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1) Determining the top of the thermocline (enter 
with Xi) 





Depth (ft) 




Depth (ft) 


1 


0 


1 


100 


2 


0 


2 


100 


3 


see X 2 


3 


100 


4 


0 


4 


100 


5 


30 


5 


150 


6 


30 


6 


150 


7 


50 


7 


200 


8 


see X 2 


8 


600 



2) Determining the bottom of the 
(enter with X 2 ) 

X 2 Depth (ft) 



1 

2 

3 

4 

5 

6 

7 

8 



300 
200 
150 
300 
300 
300 
300 
see Xi 



1 

2 

3 

4 

5 

6 

7 

8 



thermo line 



Depth (ft) 
100 
50 
30 
100 
50 
100 
100 
600 



Note that when checking for the top of the thermocline that 
when = 3 or 8, this implies looking at the second digit. 
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Likewise, when determining the bottom of the thermocline 
and X 2 = 8, then look at the first digit. 
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TABLE I 



Average Slopes Pooled Over All Possible Data 



Slope 



Source/Receiver 
Depth (ft.) 


Cruise 

No. 


di/d 

- 5 Kyds 


* 

■2 

^5 Kyds 


d2/di* 

- 5 Kyds ^5 Kyds 


50/200 


5 




.003113 


.06622 .001893 




7 






.001978 




8 


.003122 


.001446 






9 


.004476 


.001040 






10 


.003580 


.001894 






11 


.010430 


.000683 




30/150 


9 


.004950 


.000759 


.000358 




10 


.004975 


.001400 


.004202 .001356 


20/100 


5 


.028020 


.002022 


.020260 .003018 




7 




.0019581 


.002236 




8 


.005097 


.001348 


.004221 .015560 




9 


.009709 


.001062 


.008700 .001550 




11 


.009000 


.001000 


.010100 .001000 


50/400 


8 


.003831 


.001400 


.003430 .001358 




9 


.004455 


.001191 






10 


.003607 


.001924 






11 


.011110 


.000997 




100/400 


5 




.000599 






7 




.000930 






8 


.004718 


.001453 


.001964 .001313 




9 


.006369 


.001229 






10 


.003980 


.001434 






11 


.016850 


.001105 





Soxirce/receiver 



depths 
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Summary of Kolmogorov-Smirnov Goodness-of-Fit Test 
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TABLE IV. 



The Effect of Mixed Layer on Var (y ( t) ) 



Source/Receiver 
Depth (ft.) 


•k 

Above 


Crossing 


20/100 


14.8 


28.4 


20/500 


57.2 


58.7 


30/150 


19.2 


- 


50/100 


20.5 


24.9 


50/250 


5.6 


25.7 


50/400 


18.6 


28.5 


50/500 


14.6 


34.3 


100/400 


11.7 


25.4 


100/500 


12,0 


41.5 



Both source and receiver in the mixed layer 

"fc "fc 

At least source or receiver below the mixed layer, i.e., 
crossing the thermocline 
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TABLE V. 



The Gross Effect of Geographic Location 
on Var ( A db loss) - 8 Kc 



Source/Receiver 
Depth (ft.) 


Latitude 

(degs) 


Longitude 

(degs) 


Crossing 


Above*** 


20/100 


All 


All 


28.30(192) * 


14.82 (101) 




50-80N 


All 


23.12(111) 


2.60(2) 




0-50N 


0-60W 


37.60(65) 


12.29(57) 




0-50N 


60-90W 


8.30(5) 


19.16(22) 


50/100 


All 


All 


24.95(216) 


20.53(167) 




50-80N 


All 


20.00(115) 


6.77 (55) 




0-50N 


0-60W 


29.96(79) 


27.56(63) 




0-50N 


60-90W 


25.95(11) 


29.10(29) 


50/250 


All 


All 


25.71(151) 


7.53(48) 




50-80N 


All 


23.54(94) 


4.86(39) 




0-50N 


0-60W 


20.52(35) 


24.61(6) 




0-50N 


60-90W 


43.26(22) 


7.89(3) 


100/500 


All 


All 


41.58(116) 


11.97 (40) 




50-80N 


All 


39.39(87) 


12.51(34) 




0-50N 


0-60W 








0-50N 


60-90W 


49.87 (28) 


.43(3) 



*( ) - Indicates sample size. 

**At least source or receiver below mixed layer. 
***Both source and receiver in the mixed layer. 
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50-80N All 0-30 8 A 2.3(7) 7.4(17) 11.9(7) 



20/100 0-50N 0-60W 0-30 8 C 34.2(7) 44.6(21) 9.1(3) 79.7(2) 27.2(2) 

0-50N 0-60W 0-30 8 A 3.0(7) 13.7(33) 14.0(16) 
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50-80N All 0-30 8 A 2.6(2) 

50-80N All 0-30 25 C 113.8(39)117.3(30)100.8(2) 

50-80N All 0-30 25 A 54.5(2) 
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50/100 50-80N All 0-30 8 C 20.6(78) 25.5(36) 

50-80N All 0-30 8 A 8.4(41) 1.5(10) 2.5(3) 

50-80N All 0-30 25 C 31.6(44) 23.2(31) 

50-80N All 0-30 25 A 11.4(30) 2.9(8) 3.7(3) 
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1330 72/72 72/61 10 35 D 5385 

1330 83/87 - /87 20 22 D 2850 








I 




Fie 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 



TABLE X 



Columns 



1-2 

3-4 

5-6 

7-8 

9 

10 

11-12 

13-14 

15-16 

17-18 

19-21 

22 

23-25 

26-28 

29-31 

32-34 

35-37 

38-40 

41 

42-43 

44 

45-46 

47-48 

49-50 

51 

52-53 

54-55 

56 

57-60 



Acoustic Data Records 



Format 



Description 



I 2 
I 2 
I 4 

A 2 

I 2 
I 2 
I 2 
I 2 
I 3 
I 1 
I 3 
I 3 
I 3 
I 3 
I 3 
I 3 
I 1 
I 2 
I 1 
I 2 
I 2 
I 2 
I 1 
I 2 
I 2 
I 1 
I 4 



Cruise Number 
Station Number 
Hours (GCT) 

Minutes (GCT) 

Sign of Time Zone 

Time Zone 

Day 

Month 

Year 

Range (thousands of yards) 
Range (hundreds of yards) 
Receiving Ship 
South depth (ft.) 

Receiver depth (ft.) 
Propogation loss (db) 2.2 KC 
Propogation loss (db) 8 Kc 
Propogation loss (db) 16 Kc 
Propogation loss (db) 25 Kc 
Sea State ) 

BT Pattern Code ) 

Sea State ) 

BT Pattern Code ) 

Degrees (latitude) 

Minutes (latitude) 

N or S 

Degrees (longitude) 

Minutes (longitude) 

E or W 

Water depth (fathoms) 



Ship 1 
Ship 2 
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TABLE XI. 



Input Array Description - NK(I) 



Description 



Cruise number. Zero implies all cruises, otherwise 
specify cruise number desired. 

Station number. Zero specifies all stations, 
otherwise specify station number desired. 

Time (GCT) « Zero implies all times, otherwise 
specify lower bound on time interval desired. 

Time (GCT) . Same as NK(3) except specify upper bound 
on desired time interval. 

Time zone. Zero implies all time zones, otherwise 
specify desired time zone. If negative time zone, 
fill in first with negative sign, otherwise leave 
first character blank. 

Montho Zero implies all months, otherwise fill in 
lower bound no month interval desired. 

Month. Same as NK(6) except specify upper bound on 
month interval desired» 

Year* Zero implies all years, otherwise specify 
year desired. 

Range. Zero implies all ranges, otherwise specify 
lower bound on range interval. 

Range, Same as NK(9) except specify upper bound. 

Receiving ship. Zero implies both ships, 1 implies 
Ship NOo 1, 2 implies Ship No. 2. 

Source depth. Specify source depth desired. 

Receiver depth. Specify receiver depth desired. 

Latitude. Zero implies all latitudes. This is a 
packed word which expresses latitude in degrees 
and thousandths of degrees. For example: 40.817 

degrees will be put in the array as 40817, and 



TABLE XI. (Cont'd) 



14 31,7 as 31700. This will always be the lower 
(cont'd) limito If working with a n/n or S/S situation, 

enter the smaller of the two latitudes. If 
working with a n/S situation, enter the north 
latitude. 

15 N or S. This will always be an N for the n/S 
situation. 



16 



17 

18 



19 

20 



21 

22 

23 

24 

25 

26 



Latitude. Zero implies all latitudes. This is 
a packed word (see NK(14) ) . This will always be 
the upper bound on the latitude interval. If 
working with a N/N or S/S situation, enter the 
larger of the two latitude limits. If a N/S 
situation, enter the south latitude. 

N or So This will always be S for n/S situation. 

Longitude. Zero implies all longitudes. This 
is a packed word (see NK(14) ) . This will always 
be the lower limito If working with E/E or W/W 
situation, enter the smaller of the two longi- 
tudeso If E/VJ situation, enter the east 
longitude o 

E or W. This will always be E for E/W situation. 

Longitude. Zero implies all longitudeso This is 
a packed word (see NK(14) ) . This is always the 
upper limit on longitude. Enter the larger of 
the two longitudes here if a E/E or w/W situation. 
Enter west longitude for E/W situationo 

E or Wo If e/w situation, enter W here. 

Water depth. -1 implies all depths, otherwise 
specify lower limit. 

Water depth. Same as NK(22) except specify upper 
limit. 

Transmitter frequency. 13, 14, 15 and 16 imply 
2.2, 8, 16, and 25 kc respectively. 

Number of stations in cruise 5. 

Range differential between pairs of observations. 
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27 

28 

29 

30 

31 

32 

33 

34 

35 



TABLE XI. (Cont'd) 



Number of stations cruise 7. 

Niimber of stations cruise 8. 

Number of stations cruise 9. 

Niamber of stations cruise 10. 

Number of stations cruise 11. 

Number of stations cruise 12. 

Enter first cruise number to be used. 

Specify range break point. 

Number of repetitions of program desired. Must 
have a complete set of input data cards for each 
iteration. 
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TABLE XII 





Description of Arrays 


Arrays 


Description 


A 


Stores the time difference within pairs of 
observations by cruises » 


A1 


Accumulates time differences over all cruises. 


AT 


Stores time differences for pairs of observa- 
tions that have both receiver and transmitter 
neither above or below the thermoclineo 


AU 


Stores time differences for pairs of observa- 
tions that have both receiver and transmitter 
above the thermocline. 


AV 


Stores time differences for pairs of observa- 
tions with both transmitter and receiver below 
the thermoclineo 


B 


Stores the difference between adjusted trans- 
mission losses within pairs of observations 
by cruises. 


B1 


Accumulates adjusted transmission loss 
differences over all cruises. 


BT 


Same as AT for transmission loss differences. 


BU 


Same as AU for transmission loss differences. 


BV 


Same as AV for transmission loss differences. 


DBl 


* 

Decibel loss for Case I. 


DB2 


* 

Decibel loss for Case II. 


DB3 


Decibel loss for Case III. 


DB4 


* 

Decibel loss for Case IV. 


DBO 


Combines db losses for Cases III and IVo 


DB 


Combines db losses for Cases I and II. 
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TABLE XIIo (Cont'd) 



IDB 

IK 

IN 

J 

LLXl 

LLX2 

LLX3 

LLX4 

MK 

MMK 

N 

NK 

NY 

PDB 

PRNG 

RNGl 

RNG2 



Stores first and second digits for BT pattern 
codeo 

Array MK is loaded into this array one record 
at a time for computations in subroutine COMP. 

Loaded sequentially with records from the 
acoustic data tape in the main program. 

Establishes midpoints of five-minute intervals 
in subroutine SIGMA. 



Stores BT pattern 
data set in pair. 


code 


for 


Ship 


No. 1 


and 


1st 


Stores BT pattern 
data set in pair. 


code 


for 


Ship 


No. 2 


and 


1st 


Stores BT pattern 
data set in pair. 


code 


for 


Ship 


No. 1 


and 


2nd 


Stores BT pattern 
data set in pair. 


code 


for 


Ship 


No. 2 


and 


2nd 



Stores records from acoustic tape that have 
the same cruise no. 

Temporary storage for acoustic tape records 
when a different cruise number other than the 
one being considered is encountered. Will be 
placed in the MK array after computations are 
completed for the previous cruise. 

Counts entries in each five-minute interval in 
subroutine SIGMA. 

The input parameters. 

Counts entries in each 15-minute interval in 
subroutine SIGMA. 

Stores DBO for plotting at end of each cruise. 
Stores RNGO for plotting at end of each cruise. 
Stores ranges for Case I. 

Stores ranges for Case II. 
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TABLE XII. (Cont'd) 



RNG3 




Stores ranges for Case III.* 


RNG4 




Stores ranges for Case IV.* 


RNG 




"k 

Stores ranges for Cases I and II. 


RNGO 




k 

Stores ranges for Cases III and IV. 


SLOPE 


1 


k 

Stores ranges for Cases III and IV. 


SLOPE 


2 


Stores slopes for Case II* by station number. 


SLOPE 


3 


Stores slopes for Case III* by station number. 


SLOPE 


4 


k 

Stores slopes for Case IV by station number. 


X 




Stores Var (Y ( t) ) within five-minute intervals 


XX 




Stores top and bottom of the rmoc line 
information for pairs of observations. 


XRNGl 




Used in plotting RNGl. 


XRNG2 




Used in plotting RNG2. 


XRNG3 




Used in plotting RNG3» 


XRNG4 




Used in plotting RNG4. 


XRNGO 




Used in plotting RNGO. 


XRNG 




Used in plotting RNG 


XY 




Used in computing Var(Y(At)) in subroutine 
SIGMA. 


XZ 




Var(Y(At)) within each 15-minute interval 
in subroutine SIGMA. 



*Refer diagram on Page 39. 



70 
















( 







:i 

I 



1 







TABLE XIII. 



Variable 

lA 

UK 

IQ 

JJK 

JJL 

KNT 

KK 

L 

LL 

LLL 

LXl 

LX2 

LX2 

LX3 



Description of Variables 



Description 



Set equal to the station number. Used to 
check for the end of a station. 

Set to the cruise number under consideration. 
Used in finding LL. 

Counter for arrays collecting data over each 
cruise. 

May be 1 or zero. When set to zero by the 
program, it stops the iteration at the end of 
the cruise specified by NK(1) . 

A switching mechanism that allows computations 
to proceed for cruise 12 at the end of the 
acoustic data file. 

Specifies the number of iterations within each 
run of the program. 

Cruise number. 



Sets the transmitting frequency. 

Specifies the number of stations in a 
particular cruise. 



Allows for zeroing certain variables in sub- 
routine COMP-j^ at the beginning of each 
iteration. 



The BT pattern code for 
in data pair. 

The BT pattern code for 
in data pair. 

The BT pattern code for 
data pair. 

The BT pattern code for 
in data pair. 



Ship 1 and 1st set 
Ship 2 and 1st set 
Ship 2 and 1st set in 
Ship 1 and 2nd set 



71 



TABLE XIIIo (Cont'd) 



Variable Description 



LX4 


The BT pattern code for Ship 2 and 2nd set 
in data pair. 


MM 


Counter for the MK array. 


MMM 


Counter for the MMK array. 


MX 


Counter for arrays collecting data over 
all cruises. 


NC 


Set to the cruise number under consideraion 
Used to check for the end of a cruise. 


NN 


Sets the station number for use with the 
SLOPE arrays. 


SIG 


Standard deviation of Y( t) 


VAR 


Var(Y( t)) 


X 


Source depth 


Y 


Receiver depth 


Z 


Maximum allowable range differences between 
pairs of observations. 
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WRITE (6, 60) ISTA1,TIME1,TIME2tDELTIM,RNG 1»RNG2,DELRNG,AVDB1 »AVDB2 
j AV DB 

FORMAT ( '0',2X,I4,9F10.3) 
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DIMENSION X(25,l ) ,N(25) ,J(25 ) ,AA( 500) ,BB( 5C0, 1 ) ,XY( 25) ,NY(25) 
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